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ABSTRACT 



A complete reassessment of the Hubble Space Telescope (HST) observations of the transits of the extrasolar planet HD 209458 b has 
provided a transmission spectrum of the atmosphere over a wide range of wavelengths. Analysis of the Nal absorption line profile has 
already shown that the sodium abundance has to drop by at least a factor of ten above a critical altitude. Here we analyze the profile 
in the deep core of the Nal doublet line from HST and high-resolution ground-based spectra to further constrain the vertical structure 
of the HD 209458 b atmosphere. 

With a wavelength-dependent cross section that spans more than 5 orders of magnitude, we use the absorption signature of the Nal 
doublet as an atmospheric probe. The Nal transmission features are shown to sample the atmosphere of HD 209458 b over an altitude 
range of more than 6500 km, corresponding to a pressure range of 14 scale heights spanning 1 millibar to 10"' bar pressures. By 
comparing the observations with a multi-layer model in which temperature is a free parameter at the resolution of the atmospheric 
scale height, we constrain the temperature vertical profile and variations in the Na abundance in the upper part of the atmosphere of 
HD 209458 b. 

We find a rise in temperature above the drop in sodium abundance at the 3 mbar level. We also identify an isothermal atmospheric 
layer at 1 500 ± 100 K spanning almost 6 scale heights in altitude, from 10"^ to 10"'' bar. Above this layer, the temperature rises again 
to 2 500^[Q55j5 K at ~10"' bar, indicating the presence of a thermosphere. 

The resulting temperature-pressure (T-P) profile agrees with the Na condensation scenario at the 3 mbar level, with a possible 
signature of sodium ionization at higher altitudes, near the 3x10"^ bar level. Our T-P profile is found to be in good agreement with 
the profiles obtained with aeronomical models including hydrodynamic escape. 



Key words. Planetary Systems - Stars: individual: HD 209458 
Methods: observational - Methods: data analysis 



1. Introduction 

Only a few detections of extrasolar planets' atmospheric species 
are reported so far, but they have been recognized as important 
steps in our understanding of these objects. Of particular 
interest are Hubble Space Telescope (HST) observations of 
HD 209458 b during primary transit, which yielded the detec- 
tion of Nal (Charbonneau et al. 2002l l. as well as HI, OI, and 
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showing the presence of haze condensate at high alti- 
tude as previously observed in the optical using HST/ACS 
dPont et al. 20081 [Lecavelier des Etangs et al. 2008a| ). In addi- 
tion to HST, ground-based, high-resolution spectra have also 
allowed sodium detection in HD 18973 3 b (Redfield et al. 20081 1, 
confirmed sodium in HD 209458 b dSnellen et al. 20081 1. de- 
tected CO and even winds in the planetary atmosphere 
dS nellen et al. 2010b . while the other important alkali metal 
potassium has been detected in XO-2b using narrowband 



Rayleigh scattermg by H2 photometry (Sing et al. 20 10 1. 



(Lecavelier des Etangs et al. 2008b 1, along with upper limits on 
the presence of TiO and VO (" Desert et al. 20081 1 and recently 
confirmation of CII and detecti on of Silll ([Linsky et al. 2010 1 
and even possibly SilV (SchlawnetaklOlOJl. HD 189733b 
is the second transiting planet for which HST observa- 
tions yielded information on the atmospheric transmission 
specti-um. HST/NICMOS has been used to detect H2O and 
CH4 (Swain et al. 2008), though the H2O feature has been 
challenged by a more sensitive search using filter photome- 
try (|Sing et al. 20091 IDesert et al. 20091 IDesert et al. 20T0] l 



Sing et al. (2008a) revisited all the existing HST/STIS 
spectroscopic data of HD 209458 b transits gathered over the 
years. Their result is an estimate of the transit absorption depth 
(AD) as a function of wavelength from 3 000 A to 8 000 A. 
The variation of absorption depth as function of wavelength 
has been interpreted by atmospheric signatures of Na, H2 



Rayleigh scattering and TiO-VO molecules (Sing et al. 2008b 
Lecavelier des Etangs et al. 2008b , D esert et al. 2008I I. Another 



important result obtained with the analysis of this transit ab- 
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Fig. 1. The parametric T-P profile of the atmosphere of 
HD 209458 b (thick red soUd line), as adjusted over HST/STIS 
observations by Sing et al. (2008b). The three adjusted (T,P) pa- 
rameters are shown with their temperature error bars (red). The 
tripled-dot-dashed line shows the level of Na abundance drop 
at 3 mbar pressure. The thin solid line shows the 32 mbar pres- 
sure measured at the base of the sodium absorption line using the 
Rayleigh scattering. The dotted line shows the CH4/CO equilib- 
rium curve. The dashed lines show the condensation curves for 
sodium, iron, and titanium. 
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Fig. 2. The relative altitude of the absorption as observed by Sing 
et al. (2008a) close to the Nal lines in variable bandwidths (see 
text). Their photometric approach is presented as crosses with 
the coiTesponding error bars, while their spectroscopic approach 
is plotted as diamonds. The error bars coiTesponding to both 
evaluations are similar. These observations are correlated since 
extracted over variable bandwidths all centered on the Nal lines. 
Three main regions could be seen: [I] corresponding to a sharp 
decrease, [II] to a plateau and [III] to a more gentle decrease. 
Analysis of these three regions will be carried out starting from 
region [III] inward. 



sorption spectrum is the observation of a sharp drop in sodium 
abundance above the altitude coiTesponding to the absorption 
depth of 1 .485% ( |Sing et al. 20"08al|Sing et al. 2008b| l. This drop 
is clearly seen in the profile of the sodium absorption line, which 
differs from the classical "Eiffel Tower" shape found when 
there is a constant sodium abundance (e.g., Seager & Sasselov 
2000). The sodium line shows a plateau with an almost con- 
stant absorption depth around AD~ 1.485% at several hundred 
Angstroms around the line center. Closer to the line center in 
the line core, the absorption depth increases again with a shape 
more like the "Eiffel Tower" above the plateau, which shows 
that the sodium is still present above the altitude of the abun- 
dance drop. To match the width of the core of the sodium line, 
the drop in sodium abundance should be at least a factor of 10 



( |Sing etal. 20 08a Sing et al. 2008b). Moreover, if the interpre- 
tation of the AD rise at wavelengths shorter than ~ 5 000 A in 
terms of Rayleigh scattering by H2 molecules is correct, then the 
pressure at the altitude of the abundance drop is constrained to 
be about 3 mbar (Lecavelier des Etangs et al. 2008b 1. 

The fit of the HD 209458 b absorption spectrum including 
profiles of sodium and Rayleigh scattering by H2 yielded (Sing 
et al. 2008b) an estimate of the temperature-pressure (T-P) pro- 
files for the atmosphere from ~30 mbar to ~10 yubar (see also 
Madhusudhan & Seager (2009) for related issues when retriev- 
ing temperature and abundances from transit data). This result 
has been obtained using a parametric profile to fit the data, as- 
suming linear variations in temperature with altitude and two 
layers with constant slope of temperature versus altitude (Fig.[T]). 
Nonetheless, with this method, two different profiles provide an 
equally good fit to the absorption spectrum: the first profile in- 
volves Na condensation to explain the abundance drop at the 
absorption depth of 1.485%; in the second profile the tempera- 
ture is not low enough to allow the sodium condensation, and 
the decrease in sodium absorption at altitude above the plateau 
is interpreted by ionization. The two resulting T-P profiles are 



1. In the first profile the temperature decreases down to less 
than 900 K, causing a sodium abundance drop by condensa- 
tion at about 3 mbar (Fig.[T]i; 

2. The second profile (not shown in Fig. [TJ is nearly isother- 
mal at about 2200 K, in which case the abundance drop is 
explained by partial ionization of atomic sodium above the 
3 mbar level in a scenario proposed by Fortney et al. (2003). 

Both profiles do not present variations similar to model 
predictions as calculated in the literature specifically for 
HD 209458 b. The direct comparison between extracted observa- 
tional T-P profiles and model predictions will be made in Sect. 
4(Fig.0. 

To further reduce the impact of a priori hypothesis on the 
data interpretation, we propose here a new analysis of the ab- 
sorption depth profile with particular attention to the inner part 
of the Nal doublet. This new analysis is done using a non- 
parametric approach, with no a priori assumption about any ex- 
pected behaviour related to the T-P profile or to supposed Nal 
abundance variations. The hope is to be able to distinguish be- 
tween the two scenarios described above, avoiding a priori as- 
sumptions. In Sect. |2] we present the method for the data analy- 
sis. In Sect.[3]we compare the observations to the model calcula- 
tions. The results are discussed in Sect.|4]where a new descrip- 
tion of the upper atmospheric layers of HD 209458 b is proposed. 

2. The sodium absorption depth in the line core 

Sing et al. (2008a) evaluated the transit absorption depth around 
the center of the Nal lines in variable bandwidths from 4.4 to 
lOOA centered over both lines. At 12A bandwidth, the mea- 
sured absorption depth agrees with the Charbonneau et al. (2002) 
value. In the Sing et al. (2008a) approach, the absorption depth 
is evaluated relatively to a reference absorption depth in a fixed 
double spectral band (noted Ref) on both sides of the Nal lines: 
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from 5818 to 5843 A and from 5943 to 5968 A, i.e. in a fixed 
50A region. The measurements are obtained by first integrating 
the raw data over the selected bandwidth before extraction of 
the absorption depth. Because the measurements correspond to 
the total absorption depths within various bands, including the 
line center as done with filters, this absorption spectrum as a 
function of the bandwidth can thus be regarded as a spectrum 
obtained with a "photometric" approach. The resulting "photo- 
metric" spectrum is plotted in Fig.|2l 

In another approach. Sing et al. (2008a) also extracted spec- 
troscopically the absorption depth distribution over each pixel 
in the whole spectral range covered by the observations. These 
spectroscopic absorption depths can also be averaged over dif- 
ferent spectral bandwidths as in the photometric approach. We 
calculated the difference between the spectrum averaged over 
various bands centered on the Nal lines and the spectrum aver- 
aged in the reference band, Ref, defined above. The result ob- 
tained in this "spectroscopic" approach (see Figure |2| can be 
directly compared with "photometric" result. Small differences 
do show up between the two extraction methods, but they are 
well within the error bars. The main differences probably arise 
from the limb-darkening correction, which is applied at each 
pixel wavelength in the spectroscopic approach, while the cor- 
rection is a weighted average limb-darkening coiTection in the 
considered bands in the "photometric" approach. Therefore, the 
"spectroscopic" is expected to provide a better correction for the 
limb-darkening variations within the stellar lines themselves and 
in particular within the two strong Nal doublet lines. In the fol- 
lowing, we use the spectroscopic absorption depths. 

In Figure|2]the measured absorption depths have been trans- 
lated into planetary radii, i.e., altitudes above the Ref level, de- 
fined by the level of the absorption depth in the Ref spectral 
domain. To translate AD into altitudes above the Ref level, we 
used Rp=1.32Rjup jKnutson et al. 2007] l for the planetary radius 
and an absorption depth of ADuef - 1.49145% at the Ref level. 
As expected, the observed altitudes (i.e. AD) vary as a function 
of the bandwidth with a slow rise toward smaller bandwidths; 
however, three different regions could be seen in terms of AD 
variations: a region [I] (see Fig. |2]i presenting a rapid decrease 
of the AD with increasing bandwidths followed by a kind of 
plateau (region [II]), a puzzling feature discussed in Sect. |4]fol- 
lowed by a last region [III] where a slow decrease in the AD 
could be seen. These observations are not independent. In partic- 
ular, the error bars are larger than the observations' fluctuations. 
Indeed, all observations made in a given bandwidth are present 
in observation-extracted in any larger bandwidth. For that rea- 
son, we will complete our analysis of the observed behavior of 
the AD by starting from the larger bandwidths and moving in- 
ward. This approach will give us first access to the atmospheric 
information at the lower altitudes and then, by eliminating this 
information in the narrower bandwidths, move progressively up- 
wards in the atmosphere. 

3. Modeling the absorption spectrum 

The variations in the absorption depth (or variations of the corre- 
sponding altitude) along the line profile are due to variations of 
the line opacity. At each wavelength, the absorption line probes 
different altitudes in the planetary atmosphere with a resolution 
on the order of the eiTor bar on the absorption estimate, i.e. about 
0.01% in absorption depth coiTesponding to about 315 km in 
altitude. Because the absorption cross section increases toward 
wavelengths closer to the core of the lines, the absorption depth 
increases toward narrower bandwidths (Fig.|2]i. These variations 



in absorption altitude (or absorption depth) can be interpreted 
through comparison with model calculations. 

3.1. From absorption depth to physical quantities 

To interpret the measured variations in absorption with band- 
widths, we built the simplest possible model, i.e. an isother- 
mal hydrostatic uniform abundance model (IHUA). The model 
with isothermal hydrostatic atmosphere at temperature T with 
a single sodium abundance value ([Na/H], relative to the solar 
abundance), and a pressure Pj=o at a reference level, is a three- 
parameter model (r,[Na/H],Pj.=o). This model is aimed at a bet- 
ter understanding of the observed absorption depth variations in 
terms of departures from the local IHUA conditions, i.e. in terms 
of variations in temperature or sodium abundance, without any 
a priori hypothesis on the shape of the vertical temperature- 
pressure profile as needed in Sing et al. (2008a). Despite its 
simplicity, assuming hydrostatic equilibrium, this (IHUA) model 
still holds locally to interpret observational measurements be- 
cause the atmospheric scale height is on the order of the obser- 
vational accuracy, i.e. ~ 300 km for typical temperatures in the 
HD 209458 b atmosphere. The exponential decrease in the vol- 
ume density with altitude is very steep; as a result, the variations 
in the measured absorption at a given altitude are directly re- 
lated to the physical conditions (temperature, abundance) in the 
atmosphere at the corresponding altitude. 

The model absorption depth is calculated in a straightfor- 
ward manner following Lecavelier Des Etangs et al. (2008a,b). 
The total column density along the line of sight passing through 
the terminator at different altitudes is given by Fortney et al. 
(2005). The optical depth, r, in a line of sight grazing the plane- 
tary limb at an altitude z is given by 

t(A, z) ~ o-{A)n(z) ^J2nRpH (!) 

where Rpis the planetary radius, H the atmosphere scale height, 
and «(z) - n(,^o)exp(-z/H) the volume density at the altitude z 
of the main absorbent with a cross section cr{A). The scale height 
is given by the relation H = kT/fig where k is the Boltzmann 
constant, T the temperature, yU the mean mass of atmospheric 
particles taken to be 2.3 times the mass of the proton, and g the 
gravity at the planetary radius, g - MpG/Rj, with Mp equal to 
the planetary mass and G the gravitational constant. The absorp- 
tion at a given wavelength, A is calculated by finding the altitude 
at which the optical thickness is t(A, z) - Teff, where Te//=0.56 
( [Lecaveher des Etangs et al. 2008a). 

Using the equations and quantities defined above, the effec- 
tive altitude z is given by 

^ HlJ [Na/H]xP..oXcr..,(^) ^ (2) 

where cr^a/ is the Nal lines cross section. Equation |2] shows 
that the altitude (or corresponding absorption depth) depends 
on pressure and abundance, but the variation in the altitude as a 
function of wavelength does not depend upon pressure or abun- 
dance, but only upon the scale height, which is directly and 
only related to T, the temperature at the corresponding altitude. 
Therefore, the interpretation of Fig. |2] (in particular the slope 
of the absorption depth curve as a function of wavelength) pro- 
vides the temperature profile in the altitude range probed by the 
sodium line core: from 100 to 4000 km altitude. This tempera- 
ture profile can have a typical resolution close to the atmospheric 
scale height. 
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Fig. 3. The variation in the pressure as a function of wavelength 
is given according to Eq.[3](soHd Hne). Twice the solar Nal abun- 
dance is assumed. A temperature of 900K (corresponding to the 
sodium condensation temperature) is also assumed and the two 
spectral reference domains are indicated as vertical long dashed 
lines. In those regions the evaluated pressure (thick solid lines) 
is indeed on the order of 3 mbar as indicated by the horizontal 
dashed line. 



As shown by Eq. [2] the absolute values of the pressure and 
abundance impact only the absolute absorption altitudes and not 
their variations. Moreover, there is a degeneracy between these 
two quantities because only the product [Na/H] x can be 
constrained by the measurement of the absorption altitude. In 
other words, a change in the abundance can be directly com- 
pensated by a similar change in the pressure. That is the total 
quantity of gas in the line of sight. 

3.2. Pressure at the reference level 

One could assume that the pressure at the reference level, 
f;=o for our Nal study in the core of the lines, i.e. as eval- 
uated over the reference domain Ref over which the ADugf 
evaluation is made and relative to which all our differen- 
tial AD evaluations are made, is given by the H2 Rayleigh 
scattering observed at wavelengths shorter than ~5 000 A 
(Lecavelier des Etangs et a l. 2008b| l. This leads to a reference 
pressure P^=o ~3 mbar as shown in Fig.[T] This pressure corre- 
sponds to the location where the Nal abundance drops and above 
which only the core of the Nal lines do give an absorption sig- 
nature. 

To be more precise, one could alternatively use the evaluated 
Nal abundance, [Na/H] ~ 2 found in the lower layers of the at- 
mosphere by Sing et al. (2008a), to evaluate the corresponding 
pressure at the reference altitude. Indeed with only the product 
[Na/H] X P~o being defined, by using the evaluated Nal abun- 
dance in the lower layers one can deduce the pressure level sam- 
pled by the Nal lines at any spectral region. This is given by the 
following relation deduced from Eqs.[T]and|2] 



eff 



xkT 



[Na/H] X (TNaiiA) X yjlnRp/kTug) 



(3) 



The result of Eq.[3]is shown in Fig.[3]for which the assumed 
abundance is twice the solar abundance, valid up to our reference 
level. The pressure broadening in the Nal line wings is larger 




60 

Bandwidth (A) 



Fig. 4. Absorption depths calculated using IHUA models with 
r=500 K and T-l 000 K for various sodium abundances rang- 
ing from [Na/H] = 10-2, lO"-*, 10"^ 10-^ lO""' (top to bottom). 



than the natural width broadening at pressures above 29 mbar 
as it could be scaled from the corresponding values at 1 bar and 

2 000 K where the inverse lifetime related to pressure broaden- 
ing is T = 0.071cm"' X c = 2.1 X lO^s"' to be compared to the 
natural inverse lifetime which is 6.1 x lO^s ' (liro et al. 20051 1. 
Because here we are analyzing the inner parts of the Nal Unes 
(between the two reference domains, see Fig.O, we only inves- 
tigate the upper part of the atmosphere where lower pressures 
are present. We thus can ignore here the pressure broadening ef- 
fects. Consequently, cr^a/(^) in Eq. [3]was evaluated by taking 
only into account the Doppler core broadening and the natural 
wing broadening of both Nal doublet lines. 

We find that, using a temperature of 900 K and sodium abun- 
dance twice solar, absorption over the two spectral reference do- 
mains, on both sides of the Nal lines from 5 818 to 5 843 A and 
from 5 943 to 5 968 A, is due to sodium located at a pressure of 

3 mbar (Fig.|3]l. 

If a different temperature is selected for the reference level 
(as for instance 2 000 K as suggested in the second possible sce- 
nario), then the evaluated pressure at the reference level would 
be slightly higher and on the order of 6 mbar, a value still within 
the pressure domain where the pressure broadening effect could 
be ignored. With this new reference pressure, the T-P profile ob- 
tained here would still hold with only a shift in the reference 
level from 3 to 6 mbar. 

Once the reference pressure has been set, our IHUA model 
becomes a two-parameter model (r,[Na/H]) and the selection of 
a different reference pressure only slightly shifts our temperature 
evaluations along the pressure scale, with the sodium abundance 
modified by the same ratio. 



3.3. Temperature profile 

In the framework of our study, we consider relative absorption 
depths, calculated by the difference in absorption depths be- 
tween two spectral domains: the first one for a given bandwidth 
centered on the sodium lines, the second one in the constant Ref 
domain. If the [Na/H] abundance varies, the absolute absorp- 
tion depth of all spectral domains changes by nearly the same 
amount (see Eq.|2|. As a result, the relative absorption depth in 
the sodium line does not depend on the sodium abundance as 
shown in Figure H) Indeed a change in the abundance produces 
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about the same shift in the absorption depth at all wavelengths. 
Therefore, the difference between the absorption depth in a given 
bandwidth and the absorption depth in the reference wavelength 
domain (AD- AD Rgf) barely depends on the sodium abundance. 
On the contrary, the slope of the models clearly changes with 
temperature. The sodium abundance thus cannot be estimated 
by the transmission spectrum of relative absorption depth. 

However, the variations in the relative absorption depth does 
allow estimating the local temperature T at the absorption alti- 
tude. In effect, the variation in the absorption altitude with wave- 
length depends on the variation in the horizontal integrated den- 
sity with altitude. For instance, if one considers two wavelengths 
Al and A2, with coiTesponding line cross section cr{Ai) > o-{A2), 
a grazing line of sight will have the same optical thickness at 
these two wavelengths if they are at altitude zi at Ai and zi 
at A2 such that the variation in partial density of the absorb- 
ing element n{z) between zi and zi compensates for the dif- 
ferent cross sections. The optical thickness at these two wave- 
lengths are equal if n(zi)o-{Ai) - n(z2)cr(A2). Therefore, the 
variation in the observed absorption altitude as a function of 
wavelength (equivalent to variation as a function of line cross 
section) reveals the variation in the partial density at the con- 
sidered altitude. At hydrostatic equilibrium, the density vari- 
ation with altitude is determined by the temperature through 
the characteristic scale height. As a consequence, the temper- 
ature constrains the slope of the absorption altitude as a func- 
tion of wavelength. The higher the temperature, the larger the 
slope. For a given absorption feature with a known cross section 
as a function of wavelength, the temperature T can be derived 
from the measurements of the absorption altitude spectrum us- 



ing ( Lecavelier des Etangs et al. 2008a| ): 



^ ^ fig I d In o- 
k \ dA 



dz(A) 
dA 



(4) 



Figure |5] shows the expected variations for fixed abundances 
of Nal over the whole atmosphere, overplotted with isothermal 
IHUA models at different temperatures. As explained above, the 
higher the temperature, the higher the scale height in the atmo- 
sphere and thus the steeper the slope. Here, the Na abundance 
corresponding to each temperature is assumed to be constant, 
while the model representations in Figure |5] are arbitrarily ver- 
tically shifted (thus keeping the same slope at each given wave- 
length) in order to match the observations above the reference 
level. 



4. Results and discussion 

4.1. The temperature profile at low altitudes (region [III]) 

Using the method described above, we start by evaluating the 
temperature at altitudes just above the reference level (altitude 
0-225 km, bandwidths 75-110 A), with the reference level 
coiTesponding to the level where the sodium abundance drops 
sharply (~3 mbar level). The temperature of the atmosphere just 
above the reference level (see Fig.|5]l is found by the slope of the 
observed variations (see Eq. |4|. In this atmospheric layer (re- 
gion [III] as shown in Fig.|2|, the absorption depth varies slowly 
as a function of the bandwidth. As a consequence, we find the 
lower temperature close to ~ 340 K. Because this temperature 
appears to be extremely low, the assumption of uniform abun- 
dance is likely not to be valid for the scale height just above 
the reference level, a region probably still affected by the Na- 
depletion mechanism (see Sect. 14.4b . From Eq.[2] the measured 




50 80 
Bondwidth (A) 

Fig. 5. Plot of the absorption altitude as a function of the wave- 
length band. Measurements of Sing et al. (2008a) are plotted: 
photometric with crosses and spectroscopic with diamonds. The 
typical error bars of all evaluations are shown only over the 
photometric evaluations to avoid confusion. The dashed line 
shows the best-adjusted IHUA (isothermal hydrostatic and uni- 
form abundance) model to the lowest altitude observations, cor- 
responding to a temperature of 340 K. For comparison the the- 
oretical IHUA models are shown at T = 140 K, 540 K (soUd 
lines), shifted in ordinate toward the lower end of the first er- 
ror bar at 105 A in order to show the constraints on the model 
slopes (not affected by the applied shift) related to the observa- 
tional error bars. This leads to an estimate of the temperature at 
the lower levels to be 340 + 200 K. Similar errors on the order 
of + 200 K are found for each observational region of about one 
scale height (see text). Finally, the 900 K model (dotted line) cor- 
responding to the sodium condensation temperature (see Fig.lTJ 
is also shown to illustrate its incompatibility with the observa- 
tions. 



altitude is determined by the quantity of [Na/H] x cr^fl/(/l). If 
the Na abundance drops with increasing altitude, for a given op- 
tical depth through the atmosphere, the Na cross section needs 
to be larger to compensate, making the transmission spectrum 
optically thick closer to the Na core. As a result, a depletion in 
sodium causes a shallower slope (or even a plateau) in the ob- 
served absorption depth versus bandwidth. In other words, if the 
temperature is measured by the slope of the absorption altitude 
as a function of wavelength assuming a constant abundance, a 
decrease in sodium abundance with altitude implies underesti- 
mation of the temperature (Eq.|4]i. The temperature of 340 K is 
therefore a lower limit of the temperature in the layer located in 
the altitude range from to 200 km above the reference level. In 
this layer, the temperature is likely higher and the abundance of 
sodium is very likely decreasing with altitude. 

Apart from the shallow slope of absorption as a function 
of the bandwidth described above, we observe two sodium 
absorption-depth plateaus at 80 A and 55 A bandwidth. This can 
also be interpreted as abundance drops at the corresponding alti- 
tudes. We however ignore these two plateau regions because they 
are entirely immersed within the error bars domain (see Fig.|5]l. 
We thus assume locally constant abundance (at least over one 
scale height) to estimate the local temperature in the different 
layers of the atmosphere. 

The temperature profile can thus be obtained step by step 
toward increasing altitude by analyzing the observed AD vari- 
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Fig. 6. Same as in Figure |5] with IHUA models matching the 
observations up to 800 km in altitude. At the lowest altitude, 
the temperature is found to be ~ 340 K (solid line). At altitude 
above ~ 200 km, we identified layers with diff'erent temperatures 
at about 550 K (dashed Une), 910 K (dotted fine), then back to 
550 K (dashed line). A few error bars are indicated to underline 
that, in the lower altitude region [III] from to 800 km above the 
reference level, the observations are on the average incompatible 
with an IHUA model slope corresponding to ~900 K, the sodium 
condensation temperature, demonstrating that condensation has 
to take place in the planetary atmosphere. 
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Fig. 7. Same as Figure|6]for narrower bandpass, i.e. coiTespond- 
ing to regions [I] and [II] at higher altitudes from 800 km to 
4000 km. 



ations in naiTower bandpasses. When starting from wavelength 
bandwidths of 110 A, working toward shorter bandwidths, and 
assuming constant abundance, the AD variations can be inter- 
preted in terms of temperature changes as a function of altitude 
with the following characteristics (see Fig.|6]l: 

i) a lower limit of 340 K in temperature is obtained at the lowest 
altitudes, 

ii) between altitudes of 200 km and 400 km, the temperature is 
found to be about 550 K, 

iii) the AD slope constraints the temperature at about 910 K in 
the layer between 400 km and 600 km altitude, 

iv) the temperature seems to drop back close to 550 K between 
600 km and 800 km altitude. 

For each layer considered, the error bars on these temperature 
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Fig. 8. Same as in Figure |6] for the narrowest bandpass (region 
[I]). Near to the line center, the cross section is the highest, there- 
fore probing the highest altitudes. The ground-based measure- 
ments obtained by Snellen et al. (2008) are plotted with crossed 
squares together with the spectroscopic evaluations (diamonds) 
of Sing et al. (2008a). The variations of the absorption depth 
(AD) between 1.5 and 3 A bandpass are consistent with a tem- 
perature from 1 500 to 4 000 K with the most plausible value on 
the order of 2 500 K. 



evaluations are of about +200 K, similar to the error estimate 
made on the lowest temperature level at 340 K as shown in 
Fig. |5] Therefore, the observed temperature variations corre- 
spond to a marginal rise in temperature with altitude. 

Most importantly, for about three scale heights from 200 to 
800 km above the altitude where the Nal abundance is seen 
to drop sharply, the atmospheric temperatures are found to be 
within the range ~ 550 - 910 K, i.e. below the sodium con- 
densation temperature (Figure[TJ. This result favors the sodium 
condensation scenario to explain the sodium abundance drop as 
described in ( Sing et al. 2008a, S ing et al. 2008b| l. Recent global 
climate models of HD 209458 b also suggest cool temperatures 
of ~ 800 K (at mbar pressures) are possible on the night-side 
of the planet and terminator ( Showma n et al. 20091 ), further sup- 
porting this sodium condensation model. 



4.2. The temperature profile at higher aititudes(regions [I] 
and [II]) 

The temperature at higher altitudes can be estimated by looking 
for the AD variations over narrowest bandpasses. There, an im- 
portant and significant leveling off (region [II]) of the AD varia- 
tion at about 800 km altitude, between 32A and 20A bandwidths, 
is seen (note in Fig.Qits position relative to the observational er- 
ror bars). This plateau in absorption depth might be interpreted 
as a local drop in temperature but this would correspond to an 
extremely low temperature of less than 100 K ! This plateau is 
more likely related to a decrease in sodium abundance which 
will be further discussed in Sect. 14.41 

Above that long leveling off, i.e. above 800 km altitude (re- 
gion [I]), the AD increases again rapidly toward shorter band- 
widths (from 15 A bandpass and below), revealing there a rise in 
temperature. Using an IHUA model fit to the AD profile, we find 
a typical temperature of about 1500 K from about 800 km up to 
4000 km altitude (Fig.|7]l. The measurement at 4000 km corre- 
sponds to the narrowest 3A bandpass in the very center core of 
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Fig. 9. Plot of the "integrated-averaged" T-P profile (diamonds) obtained in Sect. 14.31 The vertical error bars show the pressure 
domain of the corresponding temperature estimates. The triple-dot-dashed line shows the altitude of the sodium abundance drop 
measured by ^Sing et al. 2008a Sing et al. 2008b). At pressure between 10"^ and 10"^ bar, the low temperature (T < 800K) favors 
the sodium condensation scenario to explain the sodium abundance drop. Our T-P profile is consistent with the T-P profile obtained 
by Sing et al. (2008b) using a parametric fit of the same data set (dotted line) without the addition of the latest observations at higher 
altitude levels from Snellen et al. (2008). For comparison several theoretical T-P profiles are over plotted. At higher pressures 
(P>1 mbar), models of Barman et al. (2005; dash-triple-dotted line). Showman et al. (2008; long dashed line) and Guillot (2010; 
two solid lines, one for yU = at the terminator limb and the other at // = 0.1 close to the terminator limb) are shown. At lower 
pressure levels, the Yelle (2004) theoretical T-P profile is also shown for comparison (dashed line). The dot-dashed line shows the 
vertical velocity of the hydrodynamic outflow in the Yelle (2004) model (in m/s). This velocity is negligible at pressure above a few 
10"^ bar where the hydrostatic assumption is therefore justified. 



the line in HST observations. It is noteworthy that this 1 500 K 
atmospheric layer is observed to be isothermal over nearly 6 
scale heights. Over this wide altitude range, the eiTor bar on the 
estimated temperature is rather small (±100 K). There the tem- 
perature is tightly constrained in the 1 400 K-1 600 K range. 

After publication of the HST data analysis by Sing et al. 
(2008a), new measurements of the Nai absorption depth have 
been obtained by Snellen et al. (2008) using ground-based ob- 
servations. These measurements have been obtained in narrower 
bandpasses, providing larger absorption signatures, therefore 
corresponding to absorptions at higher altitudes. These new mea- 
surements are shown in Fig. [8] together with those of Sing et al. 
(2008a). Although these various measurements have been ob- 
tained using very different techniques, and despite that ground- 
based measurements are known to be extremely difficult because 
of the Earth atmospheric perturbation, the very good agreement 



between the various data sets strengthens the confidence in those 
measurements permitting us to extend our analysis in altitude up 
to more than 6 500 km. 



On the contrary, we decided not to use the recent ground- 
based observations of Langland-Shulaet al. (2009) because their 
result suffers from uncertainties in the process of continuum fit- 
ting and subtraction of the Earth atmosphere absorption. 



We find that to match the variations of the absorption depth 
toward the narrowest 1.5 A bandwidth as measured by Snellen 
et al. (2008) , the temperature at altitudes above 4 000 km needs 
to be on the order of 2 500 K and within a 1 500-4000 K range 
(Fig. [8). 
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4.3. The "integrated-averaged" T-P profile 

An important comment here before presenting the result is to 
recall that the evaluated T-P profile is extracted from quantities 
"integrated" along the cords passing through portions of the limb 
presenting gradients, e.g. from the night to day sides, as well 
as "averaged" over the planetary limb, e.g. from equatorial to 
polar regions. In a first step, such an "integrated-averaged" T-P 
profile could be compared to 1-D model calculations, although 
it is clear that later steps, when observations permit it, would be 
to compare them to more sophisticated 3-D models as illustrated 
in e.g. the Iro et al. (2005) and the recent Fortney et al. (2010) 
studies. 

From the temperature-altitude relation obtained in Sects. |4TT] 
and l4.2l and assuming hydrostatic equilibrium, we can calculate 
the temperature-pressure (T-P) profile. We start by using a ref- 
erence pressure of 3 mbar at the reference km altitude. The 
pressure is calculated step by step, starting from the zero alti- 
tude level. The result is shown in Fig.|9l 

Acknowledging that the sodium condensation occurs at tem- 
peratures below 800-900 K, the temperatures measured in the 
layer between and 800 km above the reference level are found 
to be low enough to allow condensation. This condensation 
seems to take place at pressures between 3 mbar and 10 //bar. 

As already mentioned in the introduction, at pressures above 
3 mbar, the evaluated observational profiles do not present vari- 
ations similar to the model predictions as calculated in the lit- 
erature specifically for HD 209458 b {e.g. Barman et al. 2002; 
BuiTows et al. 2003; Iro et al. 2005; Barman et al. 2005; Fortney 
et al. 2008; Showman et al. 2008; Showman et al. 2009; Guillot 
2010; Fortney et al. 2010). Three of them (Barman et al. 2005; 
Showman et al. 2008; Guillot 2010) are presented in the lower 
part of Fig. |9] It is interesting to note that all temperature pre- 
dictions at about 3 mbar are in good agreement with our evalu- 
ation simply made from the direct estimation of the scale height 
at these levels. We also point out the large variations revealed 
by the two near terminator-limb (i.e. near the transit observed 
limb) model calculations from Guillot (2010; one for // = at 
the terminator limb and the other at = 0.1 close to the termina- 
tor limb) as also shown in Fig.|9] This may indicate that indeed 
such model predictions are difficult to compare to observations 
according to their extreme sensitivity to the precise limb geom- 
etry and, as already said, to the integrated nature of the observa- 
tions along a grazing line of sight. Also, one could argue that at 
higher pressures, it is not surprising that observations (and their 
interpretation) differ from theoretical models, since for instance, 
a large heat advection could be critical in the troposphere leading 
to strong departures from equilibrium temperature profiles. 

Above this condensation layer, a progressive rise in the tem- 
perature with altitude is needed to explain the profile of the 
sodium line toward narrower bandpasses. The temperature in the 
upper atmosphere is found to reach ~ 2 500 K at the very low 
pressures of ~ 3 to ~ 90 nbar. This is likely the signature of the 
base of the thermosphere, which is the expected transition be- 
tween the cool lower atmosphere and the hot upper atmosphere 
where hydrodynamic atmospheric escape takes place. This sce- 
nario agrees with theoretical models of the atmospheric escape 
as published by, e.g., Lecavelier des Etangs et al. (2004), Yelle 
(2004) over-plotted in Fig. |9] Tian et al. (2005), Yelle (2006), 
Munoz (2007), and Murray-Clay et al. (2009). 

A summary of our T-P profile is given in Table [T] The scale 
heights H are given for each temperature layer. This shows that 
the layers have thicknesses of at least one or more scale heights. 
With a pressure change by a factor of 2.7 over one scale height 



Table 1. Evaluated atmospheric parameters 
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under hydrostatic equilibrium, Fig.|9]shows that the bottom layer 
at the lowest temperature of > 340 K has a thickness of nearly 
2 scale heights. The three following layers at 550 K and 910 K 
are over about 1 scale height, while the atmospheric layer with 
a constant temperature at ~ 1 500 K extends over nearly 6 scale 
heights. Finally the last upper layer at about 2 500 K extends 
over more than 3 scale heights. 

The calculated T-P profile is valid when assuming local hy- 
drostatic equilibrium. This assumption could be invalid because 
of the hydrodynamical escape of the upper atmosphere. Figure|9] 
shows that the theoretical escape velocities triggered by the hy- 
drodynamic flow as evaluated by Yelle (2004) are negligible be- 
low the altitude corresponding to a pressure of ~ 10"^ bar. This 
shows that the local hydrostatic equilibrium assumption is valid 
over our studied range, even at the highest altitudes where we 
estimate a temperature of about 2 500 K. 

4.4. The Na abundance profile 

The sodium abundance [Na/H] cannot be constrained in abso- 
lute terms according to the mentioned degeneracy between abun- 
dance and pressure in the absorption spectrum (see Lecavelier 
Des Etangs et al. 2008a and Eq.|2]i. However, here the degener- 
acy can be solved by assuming a reference pressure at the refer- 
ence level (see Section [3T2l i. Indeed we used a reference pressure 
of 3 mbar for the altitude corresponding to our reference level 
at 1.49145% absorption depth. Using in addition the sodium 
abundance below the reference level as evaluated by Sing et al. 
(2008b) in the lower layers of the atmosphere, [Na/H] -3x10"'' 
(solar abundance equals 1.5x10"^, Asplund et al. 2006), we can 
deduce the sodium abundance evolution in the upper layers. In 
agreement with the condensation scenario, the sodium abun- 
dance drops just above the reference level and then, following 
the upward evaluation of Sing et al. (2008b), it should be re- 
duced on the average over all the upper layers by about a factor 
10. 

Starting from the value at the reference level, the Na abun- 
dance is evaluated at higher levels using a step-by-step proce- 
dure with local IHUA models. A constant [Na/H] abundance in 
the isothermal layers obtained in Sects. I4.1l and l4~2l has to be 
present, in agreement with om IHUA hypothesis. 

A plateau in the measured absorption depth is observed at an 
altitude of about 800 km. This plateau is a kind of anomaly: with 
the increase of absorption cross section in the core of the line, 
the AD should rise toward smaller bandwidths. If the plateau 
was extending to the smallest bandwidths at the center of the 
line, this would reveal the absence of the absorber at altitudes 
above this plateau. However, here the plateau has a limited ex- 
tension from AA\ ~ 15 A to AA2 ~ 30 A bandwidth. The absorp- 
tion depth rises clearly again at bandwidths smaller than about 
15 A (see Fig.|7]i. This reveals a drop in the sodium abundance at 
the altitude of the plateau, but the sodium is still present above. 
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Fig. 10. Summary of the Na abundance estimates using the IHUA 
model. Because our approach is local we can only identify here 
the decreases in sodium abundance by a factor of about 4 at 
3 X 10"^ bar (dot-dashed line), although NLTE effects could per- 
turb that estimation (see text). The drop in abundance by about a 
factor of 10 showed at the condensation level (at 3 mbar, tripled- 
dot-dashed line) is taken from Sing et al. (2008a,b). 

thereby producing the observed absorption rise toward the line 
center. Along the plateau, the variation on the line cross section 
with wavelength, from criA^) to cr{A2), is compensated for by the 
variation in the sodium abundance, [Na/H]i and [Na/H]2 above 
and below the plateau altitude. Using Eq. |2l the altitude of the 
plateau z(/l) being constant from Ai to /I2, we find that 



aNa/H]2 



<t(A2) 
ct{A,) 



(5) 



l(r{A2)\ 




-Ao\ 


\cr{A,)] 




-Aq) 



Because here the absorption is dominated by the damping wings 
of the sodium line, we have 



(6) 



With A/l2/A/ii ~2 for the plateau at 800 km of altitude, we can 
determine that this plateau may be explained by a drop in sodium 
abundance by a factor [Na/H]2/[Na/H]i ~ 4. According to our 
approach, which is essentially local, we can conclude about this 
relative local drop, while as stated earlier, absolute values of the 
sodium abundances could not be assessed. To represent these 
absolute values, we used the result of the global approach of 
Sing et al. (2008a,b) as shown in Fig.fTOl 

Considering the sodium abundance drop evaluated by Sing 
et al. (2008a,b) at the condensation level to be on the order of 
a factor of 10, we assumed that [Na/H] 2 should be on the order 
of 0.2 solar, while the Na abundance above the altitude of the 
plateau should be [Na/H]i~ 0.05 solar. Figure [TOl shows the Na 
abundance variations obtained from the observed AD over the 
whole domain sampled in Fig.|9] 

To match the observed AD plateau, one has to assume a local 
reduction in the Na abundance by a factor of about 4, which 
could indicate of a possible physical effect in the atmosphere, 
such as the ionization of Nal into Nail as predicted by several 
models, e.g. Fortney et al. (2003). The conclusion on the drop 
in the Na abundance to interpret the AD plateau is expressed in 
terms of relative abundance, so this conclusion does not depend 
upon any assumption about pressure or absolute abundances. 

We are evaluating variations in the Na abundance in low 
volume density regions where nonlocal thermodynamic effects 



(NLTE) should take place, i.e. below 1 mbar as evaluated by 
Barman et al. (2002). There radiative processes dominate colli- 
sional ones and the Nal population levels could be significantly 
away from the assumed local thermodynamic equilibrium (LTE) 
ones. As long as such effects do not vary significantly over one 
scale height, our temperature evaluations should not be affected 
by these NLTE effects. On the contrary, such effects certainly 
should perturb abundance evaluations, and in particular, the pre- 
viously mentioned factor 4 drop of abundance at 3 x 10"^ bar 
could as well be the signature of such NLTE effects at these 
altitudes, leading probably to an underestimation of the evalu- 
ated drop (Barman et al. 2002). If the coiTesponding signature 
is indeed due to Na ionization, it is in reasonably good agree- 
ment with the Barman et al. (2002) evaluation who found that 
Na would be mostly ionized in the limb for P < 10"^ bar and 
only partially ionized (~ 3%) below this point. 



5. Conclusion 

Using the measured Nal absorption depth as a function of the 
bandwidths, we estimated temperature variations as a function 
of altitude. Our temperature profile is inconsistent with one of 
the two scenarios of Sing et al. (2008b) (scenario with constant 
temperature and ionization at low altitude). 

We found layers at temperatures clearly below the 
Nal condensation temperature showing that sodium con- 
densation has to take place at these levels, followed by a 
temperature rise at higher altitudes up to about 2500 K. 
This rise is likely the signature of the thermosphere 
needed to link the cool lower to the hot upper atmo- 
sphere, as suggested by the observations of atmospheric 
escape 



Vidal-Madjar et al. 2004 



(|Vidal-Madjar et al. 2003 

Linsky et al. 2010|l, observations of high temperatures 
(Ballester et al. 2007]|, and models of the escape mechanism 
(Lecavelier des Etangs et al. 2004[ lYelle 20041 ITian et al. 20051 
lYelle 2006ilMunoz 2007i|Murray-Clay et al. 2009| ). 

The results of our analysis favor the sodium condensation 
scenario to explain the deficiency of sodium observed in the 
line core, and the possible presence of an Na ionization layer 
just below the base of the thermosphere. Our temperature profile 
presents the following patterns: 

- a rise from 340 K to 1 500 K over -800 km (about 5 scale 
heights) in the pressure range of 3 x 10^^-3 x 10"^ bar; 

- a thick layer with a nearly constant temperature at about 
1 500 K. This layer has a thickness of about 3 200 km cor- 
responding to nearly 6 scale heights for a pressure range of 
3 X 10-^-10-^ bar; 

- a further temperature rise is observed at higher altitudes, 
up to 2 500 K at the highest altitudes about 6 500 km above the 
reference level, corresponding to a pressure of about 3 x 10"^ bar. 
This high temperature is likely related to the bottom part of ther- 
mosphere whose presence is predicted by theoretical models of 
the atmosphere. 

- an Nal abundance drop by a factor of 4. This is possibly due 
to ionization at the 3 x 10"^ bar level, the strength of the observed 
drop depending upon the importance of the NLTE effects at these 
altitudes. 
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